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Abstract 

We calculate thermal corrections to the non-linear QED effective action for 
low-energy photon interactions in a background electromagnetic field. The 
high-temperature expansion shows that at T ^> m the vacuum contribution is 
exactly cancelled to all orders in the external field except for a non-trivial two- 
point function contribution. The high-temperature expansion derived reveals a 
remarkable cancellation of infrared sensitive contributions. As a result photon- 
splitting in the presence of a magnetic field is suppressed in the presence of 
an electron-positron QED-plasma at very high temperatures. In a cold and 
dense plasma a similar suppression takes place. At the same time Compton 
scattering dominates for weak fields and the suppression is rarely important 
in physical situations. 
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1 Introduction 



Recently the non-linear effects induced by virtual electron-positron pairs in quantum 
electrodynamics (QED) have been discussed in great detail in the literature ([I], |3| and 
references cited therein) confirming earlier calculations concerning in particular photon 
splitting in external magnetic fields ||, H. It is of great interest to notice that such non- 
linear QED effects may lead to observable physical effects in e.g. the gamma-ray burst 
spectra (for reviews see e.g. H 0]) in particular with regard to the so called soft gamma- 
ray repeaters ||, [|. Recently, photon splitting processes in the magnetosphere of 7-ray 
pulsars has also been discussed |IJ , where it has been argued that such processes can be 
comparable to pair-production processes if the magnetic field B is sufficiently high, i.e. 
B > B c , where B c si4x 10 13 Gauss is the critical field in QED (i.e. eB c = m 2 ). 

In vacuum the presence of virtual electron-positron pairs leads to the well-known 



Euler-Kockeljl l |-Heisenberg[|l2|]-Weisskopf [ 13 1 — S chwinger ||T4|1 (EKHWS) effective action, 
which was used in the classical paper by Adler on the subject of computing photon- 
splitting amplitudes and absorption coefficients in magnetic fields || . Astrophysical mod- 
els of neutron stars, which have been used in explanations for the attenuation of 7-rays, 
suggest that photon-splitting processes in the presence of strong magnetic fields do not 
necessarily take place in an environment which may be approximated by a vacuum. In- 
stead finite temperature T and/or chemical potential fi may have important effects which 
motivates a study of photon splitting at finite temperature (see e.g. Jl5|, [L6|, |l7], [HI)- We 
are considering the limit of weak fields in order to see how large the effects can be and it 
turns out that for an electrically neutral e + e~-plasma in a weak field the Compton scat- 
tering is dominating whenever the thermally induces splitting amplitude is appreciable. 
There is a possibility that the thermal splitting rate becomes larger than the Compton rate 
for finite chemical potential and very low temperature, but our calculational technique 
breaks down in that limit. 

It has been noticed in the literature that if T > m all the non-linear terms in the 
effective EKHWS action are cancelled by thermal corrections and the remaining effective 
action becomes quadratic in the electro-magnetic field strengths with a non-trivial depen- 
dence in T, ji |L9|, |2(J ^1], |22j. Below we reconsider this cancellation explicitly including 
0(m/T) corrections. As a function of the photon energy u the results obtained in the 
present paper strictly only apply to the situation where w/m< 1. In the vacuum sector 
Adler || has, however, calculated the exact u/m dependence and verified that this de- 
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pendence on the photon splitting processes is rather weak apart from phase-space factors. 
The u/m- dependence can therefore be extrapolated to oj/m ~ 1, at least for magnetic 
fields such that B/ B c < 1. We expect a similar weak a;/m-dependence also in the presence 
of a thermal environment even though this has not yet been verified explicitly. Since we 
are mostly interested in the temperature and magnetic field dependence we will simply 
put to = m. 

The paper is organised as follows. In Section |2| we recall the form of the QED effective 
action including one-loop thermal corrections due to the presence of a fermion-antifermion 
heat-bath. In Section [| we perform a weak field expansion of the effective action and a high 
temperature expansion is considered in Section |j. In Section [| a similar low temperature 
expansion is given. A non-covariant contribution to the effective action, fourth order in the 
external fields, together with sixth order vacuum and thermal corrections, are considered 
in Section ^ with regard to their effect on photon splitting processes. Final comments are 
given in Section [7|. 



2 Effective Action 



The thermal one-loop non-linear QED effective Lagrangian £^(E,B) for slowly varying 
electric and magnetic fields has been calculated |20| with the result 
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The field-independent contribution of Eq. ( [2.1| ) can be obtained by making use of dimen- 
sional regularisation techniques as described in JH]]. The function h(s) = (eFcotheFs)oo, 
where F^ u = d^A u — d u A^, was calculated in |2(J and is given by 
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In terms of the field-strength we use the notation 
E 2 + B 2 
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where the Lorentz-invariant quantities T and Q are given by 

,2 _ b 2 B 2_ E 2 



a 



Q = ab = E • B 



(2.4) 



in terms of the magnetic and electric fields B and E (B = |B|, E = |E|). In |2(J it was 
argued that the fermion equilibrium distribution function, _/>(po, eA ), in the case of an 
external electromagnetic field, should be chosen in the following form 



f F (p ;eA 



6(p - eA ) 9(-p + eA ) 



*P(po~IJ) _|_ 1 e /3(-Po+M) _|_ 1 ' 
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where j3 = 1/T is the inverse temperature and [i the chemical potential related to the 
conserved charge of the system. In this notation the renormalised low energy EKHWS 
vacuum effective action for non-linear electromagnetic fields in QED has the well-known 
form (see e.g. |23| ): 



£ vac (E,B) 



» ds 

87I" 2 Jo S 3 



s 2 e 2 



s e afecoth(esa) cot(esfe) — 1 — (a 



(2.6) 



3 Weak Field Expansion 



The total effective action for a slowly varying background field is given by 

£ eff (E, B) = —T + £ vac (E, B) + £&"(E, B) . 



(3.1) 



The vacuum contribution £ vac (E, B) has a straightforward expansion in terms of the 
invariants T and Q. We are now interested in a similar weak- field expansion of the 
thermal contribution /^'^(E, B). Care must then be taken since, as we will shortly see, 
one easily encounters severe infrared problems. The basic idea is to collect everything in 
the integrand in Eq. ( |2.1|) which depends on the fields into one function multiplying the 
zero- field part and then expand it formally in E and B. For a strictly degenerate plasma 
£^ ,Al (E,B) is not analytic in the fields but shows de Haas-van Alphen oscillations as 
discussed in f2~I|| . However, at finite temperature and weak enough fields these oscillations 
average out over the Fermi surface and a weak field expansion becomes meaningful. In 
this situation we can neglect the ie in the combination h(s) — ie which is only needed for 
the non-analytic structure. In particular we need the following expansion 

oo 

= Y,ak~2(uj;H;F;g)s k , (3.2) 
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where u = po — eA turns out to be a convenient variable of integration. The coefficients 
ctk(uj; H; T\ Q) are rather lengthy and we have most easily evaluated them using a symbolic 
computer program. After performing the s-integrals in Eq. ( |2.1| ), using the expansion in 



Eq. ( |3.2| ), we find that 



£^(E, B) = -—^ 2 J ^ duf F (u)Im £ «*H (-V - m 2 ) + e) 2 r(* + -)e** 

(3.3) 

where we in ak{oS) have suppressed the dependence of the variables 7i, T and Q, and 
where we have defined /f(<^) = fp{u + eA ;eA ). When making use of the variable 
u, the thermal distribution function in Eq. ( |2.5| ) is naturally expressed in terms of the 
effective chemical potential fi e g = fi — cAq. If we in Eq. ( |3.3| ) took the e — > limit at 
this point we would get serious artificial infrared problems with the cu-integral. We avoid 
these divergences by first performing a number of partial integrations using 

(^-m=) +£ )-"'' = - ;T |^ y ^(- J (^- m 2 ) + E )- ,/2 , (3.4) 

where is a derivative with respect to u 2 . From now on we shall drop the k = —2 
term since it only gives the field independent part (i.e. the free energy in the absence 
of the external fields) which comes out in a standard way \2~I\. For the rest of Eq. 



denoted by AC 13 '^, we rewrite D^ 2 as derivatives with respect to id and integrate them 
all by parts. This is not a problem since the thermal distribution function is zero at the 
boundaries. We obtain the result (notice that ak=~i(uj) = 0): 



A£P>" = -J- [°° dw ^ f m ^ fy ( Al)\a fc (,)/ F M) , (3.5) 
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where we have used Im (— i{oj 2 — m 2 ) + e) _1 / 2 e^ — > Q(uj 2 — m 2 )(uj 2 — m 2 )" 1 / 2 . The 
expansion in Eq. (|3.5| ) is only valid if /f(w) is smooth enough. Otherwise, a non-analytic 
structure like the de Haas-van Alphen oscillation can occur. We can then write 

A£f>>» = - JL r d J { ^~ m2) f] e 2n A^C^(u) , (3.6) 
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where A^C^ (u) is of the order 2n in the electromagnetic field-strengths E and/or B. 
The lowest order terms are given through 

A . ¥fH= W, + 4"H^ (37) 

3 6 
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and 



60 a; 4 



where we have used Ji — T = E 2 and the notation fpioj) = d n f F {u) / duo 11 . 



4 High Temperature Expansion 

In the very high temperature limit, i.e. T ^> m, |// e ff|, we see that /f(u;) — > | and 
therefore only a few terms in the expansion of A( 2n )£^(u;) survive. The quadratic terms 
go like log(T/m) at high temperature and the higher order corrections go to constants. 
Then the cu-integral can be performed for the higher corrections using 

± / ( " 2 -'"V * ^> (4 1) 

\/uj 2 — m 2 Xm^F^^) 

where k > . In this limit the effective Lagrangian then becomes 

e 2 T e 2 TP 2 

AC^ = -—JE 2 -B 2 )H-) + e -^ 



360vr 2 m 4 

" 1260^ [(jg2 " ( 2(jg2 " + 13(E " B)2 ) ] + ° {F8) ' (4 ' 2) 
where F generically stands for an external field. It is interesting to note that the higher 
order corrections are exactly equal, but with opposite sign, to the vacuum corrections 
obtained from Eq. 



n , c(vn] (a 2 -b 2 ) 2 + 7(ab) 2 (b 2 -a 2 )(2(b 2 -ar + 13(ab) 2 ) 

C (E ' B) = 360^ + 1260^ + <D{F) ■ (4 ' 3) 



As a matter of fact, by inspection of Eq. (|2.1|) and Eq. ( |2.6|) , this cancellation is true 
to all orders in the external fields except for quadratic terms exhibited in Eq. (|4.2| ) |20| . 
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In passing, we also notice that at zero temperature and for a large chemical potential, 
l/ieffl S> m, the thermal non-linear effective action CP ,tJ, (E,B) also cancels the vacuum 
effective action C vac (E,B). 

At high temperature the vacuum and thermal effects thus cancel and we have to go to 
higher order in m/T to find the leading contribution. If we first subtract the f F {oj) = | 
part and the scale out a factor of T~ 2n ~ 2 in the remaining 0(F 2n ) terms, we find that 
the integrand is infrared finite as m/T — > 0. This is far from obvious and comes from a 
remarkable cancellation of different infrared sensitive terms. To get the first sub-leading 
term we thus take the limit m/T —>■ and then we are left with a dimensionless integral. 

To be more explicit let us take /x cS = and write the subtracted 0(F 4 ) Lagrangian 

as 



A£ (4 )= i rdp 



2ir 2 Jo uj 



- 4 (f F H - ~) + ^D u f F {u>) 
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+ ^l D 2 f F (uj) + ^Dlf F {u) + CDtMuj) 
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(4.4) 



where now f F {uj) = (e^ + l)" 1 and p = Vuj 2 - m 2 . The full expression is perfectly finite, 
due to cancellation mechanism alluded to above, but each term diverges at (3uj = 0. To 
be able to treat the terms separately we use a zeta-function regularisation technique by 
multiplying the integrand in Eq. (|4.4j) by uj u , considering v large enough first, and then 
take the limit v — > at the end of the calculation. Let us therefore define the function 

'<*■-> s j£ wfw • (4 ' 5) 

and notice that 

I(z,n) = I(z,n-l) + ^-I{z + l,n-l) , 

n — 1 

v ' ; 2 Z sin(z7r/2) v ; 

For positive z the function I(z, n) has been defined in terms of an analytical continuation 
using the Riemann (^-function. The integrals in Eq. (4.4) can then be performed by partial 



integrations, where the boundary terms vanish for large enough u, and thereby reducing 
the result to 

A£W = ^|i[ 54 + 6 ^ 2 - 12 ^ + 9 ( E - B ) 2 ] • ( 4 - 7 ) 
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The 0(F e ) correction can be treated in a similar way and we find 



511C(9) 



147456vr 10 T 8 
+279E 2 (E • B) 2 



[2B b + A2B 4 E 2 + 3QB 2 E* - 220E 



7B 2 (E-B) 2 



(4.8) 



5 LOW TEMPERATURE EXPANSION 



In many applications, such as photon splitting (see Section ||), it is more interesting to 
have a low temperature expansion. At high temperature Compton scattering would any- 
way dominate over pure splitting (at least for low eB/m 2 ) and in several astrophysical 
situations the temperature is relatively low compared to the electron mass. On gen- 
eral grounds we expect low temperature effects to be exponentially small since they are 
Boltzmann suppressed, but as we shall find there is a region where the effects still are 
appreciable. The formula we need for the expansion is for fi e g = 
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We find that, although all terms go to zero exponentially, the ones with more derivatives 
with respect to u have higher inverse powers of T which makes them relatively large for 
T/m < 1, where also the approximations in Eq. (|5.1| ) are numerically good. Using these 
terms the effective action to 0(F 4 ) becomes 
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(5.2) 



6 Weak field photon splitting 



In the classic paper by Adler || the QED photon splitting process in a background 
magnetic field is treated in great detail. One important observation in that paper is that 



s 



the 0(F 4 )-term does not contribute to the splitting amplitude in vacuum. This result 
relies on the covariant form of the effective action, i.e. that only T and Q occurs in 
the explicit expression, and that all the photons are collinear due to kinematics. Non- 
linear effects on the dispersion relation either make it possible to have a non-zero opening 
angle between the split photons or prohibits the splitting altogether, depending on the 
polarisation states. At finite temperature the situation is, however, different since now 
the manifest Lorentz invariance is broken in the effective action by the dependence on 7i. 
The consequence is that photon splitting is allowed already from the 0(F A )-term even for 
collinear photons. 

Thermal corrections also alter the dispersion relation of the photons which changes the 
kinematic conditions for splitting, and then not all combinations of photon polarisation 
are possible. These effects were studied in | on a general basis and it was found that 
the only allowed process is ||— >J_i + _l_2, provided the electron density is not too large 
(n e < 10 19 cm~ 3 ). Thermal corrections to the dispersion relation have several origins. 
First, there is the plasma mass which is independent of the background field. The static 
version, the Debye mass, can found from Eq. ( ^.1[ ) by noting that the O(F ) term depends 
on A through the thermal distribution function. Expanding O(A^) gives the Debye mass 
||20|| . Secondly, higher order terms like the JF?i-term change the dispersion relation in 
a background field. Details of the dispersion relation are complicated by the fact that 
polarisation tensor is non-analytic in the energy and momentum. We shall concentrate 
on situations where the photon frequency u is smaller than the electron mass so that 
our effective action for slowly varying fields is valid. At the same time we take the 
temperature and chemical potential to be low so that correction to the dispersion relation 
is small and the photon almost light like. In this way we do not have any large correction 
to phase-space integrals which would arise when the thermal mass is important. To be 
more concrete we first consider an e + e~-plasma with /i e fj- = and T/m sufficiently small. 
Then we require that 



which in actual numbers means T/m <C 4.5. It may at first seem inconsistent to assume 
the effect on the dispersion relation to be small and still argue that plasma effects on 
the splitting amplitude can be important. However, corrections to the dispersion relation 
should be compared with the photon energy and momentum, while thermal splitting 
amplitudes should be compared with other splitting or scattering processes. Moreover, 




(6.1) 
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higher order terms in the field strength are more IR-sensitive at low T and therefore 
relatively larger. 

The new ingredient in a background plasma is the FTi-term and the question is 
whether there is any regime in which the thermal splitting is important if the temperature 
is not much higher than the electron rest mass. At low temperature the vacuum effects 
dominate and at high temperature Compton scattering is overwhelmingly large. 

The various processes to be considered depend differently on the photon energy uj and 
the scattering angle 9 between the direction of the propagation of the initial photon and 
the external magnetic field, but we shall choose for reasons of simplicity sin# = 1 and 
uj = m in the final numerical estimates. This choice of u is in the upper limit of the 
validity since our approximation assumes a low photon frequency (almost constant fields) 
and it should be remembered that the splitting probability goes like (oj/m) 5 (see below), 
and thus decreases rapidly for smaller uj. 

In the low photon frequency limit the amplitude for the ||— >J_i + J_2 photon splitting 
process is given by 

M[|H± 1 + J 4 ]=^{((fcxgH) i A + ^l»] 




e,b: 



3 ^^3 , 

The photon absorption coefficient is then given by 



E = 



k(|H± 1 + ± 2 ) = — — 2 I ' duj x I ' duj 2 S(uj-uj 1 -uj 2 )\M[\\^± 1 + ± 2 ]\ 2 . (6.3) 
OZ1T(jJ z Jo Jo 

The appropriate vacuum 0(F 6 ) splitting probability for the process ||— + _L 2 can be 
obtained from Eq. ( [4.31) : 

feB\ 6 /u\ 5 Rn 169a 3 , . 

i — ■ ' — 1 sin 6 ^— — — . (6.4) 



(6.2) 
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The only difference from the well-known result of is that we use units where a = e 2 /47r 
and = V^F^. 

Next, we need the splitting probability from the thermal 0(F A ) term. The non- 
covariant part of C 13 ^ to that order, which is the only one contributing to photon splitting, 
is obtained from Eq. (|3.8| ) and is given by 



A „ fl „ 1 B 4 too e(u 2 -m 2 ) d ff {2) (uj)\ , x 

AC^ = - — — du ) > —lLJU.\ . 6 .5 

2n 2 360 J-oo y/uj 2 — m 2 du \ uj J 
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As in the case of photon-splitting processes in the vacuum, we derive the splitting rate 
for the process ||— >J_i + _L 2 using of Eq. ( |6.5| ) 



«J ^.) 2 («)\^e^Ur~( f -^)) 2 • (6.6) 



m \m* J \mj 303757T 2 y jo uj du \ uj 

To be consistent we also compute the thermal 0(F 6 ) splitting rate Kg in a similar way. 
Finally it should be noticed that the amplitudes for all these processes should be added 
coherently since the final states are the same and the external field is constant and thus 
coherent. The total splitting rate K tot shows the characteristic decrease for high T when 
the thermal contribution interfere destructively with the vacuum part. The three contri- 
butions and the coherent sum are compared in Fig. [I]. (The physical attenuation length 
in a QED plasma can be obtained using m^ 1 = 3.86xl0 _11 cm.) 

In addition to splitting the photon can scatter directly with the plasma which turns 
out to be the dominant process for a large parameter range. We estimate the scatter- 
ing rate using the total Compton cross-section which for unpolarised photons is ac = 
87ra 2 f(uj/m)/3m 2 , where f{uj/m) is a slowly varying function such that f{uj/m) ~ 1/2 if 
uj — m (see e.g. p3[). In order to get the absorption coefficient we multiply oq with 
the appropriate density of electrons and positrons 

™ = ?£i 2 r dp 1? . (6 .7) 

m m 7r 2 Jo (eft* + 1) 

It is clear from Fig. [1] that the only region where the thermal splitting can dominate 
over the vacuum one it is already far below the Compton scattering rate and will thus 
not be very important. 

As a second example we shall consider the, for astrophysics, more realistic case of an 
electron-proton plasma. Since the electron is much lighter than the proton it dominates 
the correction to the splitting amplitude so we shall only take into account the electron 
part. For a non-degenerate plasma at relatively low temperature (T <C m e ) the effective 
chemical potential is close to the electron mass (larger /j, e s gives a degenerate plasma 
and smaller /i e g- corresponds to exponentially small electron density). In order to see 
the main features we can therefore put fi c g = m e . In Fig. |2| we can see that now the 
Compton damping rate rapidly decreases for low T while the splitting rates, which are 
more IR sensitive, actually increase. The fact that the thermal splitting rate becomes 
large shall in this case be taken as an indication of that we approach the limit where 
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0.0 0.5 1.0 1.5 2.0 



T/m 

Figure 1: The vacuum and thermal photon splitting absorption coefficients ACg ac , 
/C4, Kq , and the Compton absorption ko (upper solid curve) are shown as a function 
of temperature in units of m for eB = 0.2m?, uj = m and sin# = 1. The lower solid 
curve corresponds to the coherent sum of amplitudes up to 0(F 6 ). By including 
thermal corrections of 0(F e ), the splitting amplitude goes to zero at sufficiently 
large T/m as discussed in Section ||. 



the weak field expansion is no longer valid. We know that for T = and // e g > m e we 
have de-Haas-van-Alphen oscillations which are non-analytic in B and that is reflected 
in divergent coefficients in the weak field expansion. Also, k§ diverges more rapidly than 
k± which shows the breakdown of the perturbative expansion f2l| . On the other hand 
Fig. ^| also tells us that there might be a region for realistic plasmas where the thermal 
splitting actually is the dominant process, but to settle the question it would be necessary 
to perform the calculation using the exact Landau levels without expansion in powers of 
B. 
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0.0 0.2 0.4 0.6 0.8 1.0 



T/m 

Figure 2: The same symbols and values of parameters as in Fig. [l] except that 
here the chemical potential is /i e fj = m. 

7 Final comments 

A detailed analysis of the astrophysical consequences of thermal photon splitting is out- 
side the scope of this paper but we have found that in many generic cases this process 
is subdominant and not likely to be important. We have shown in Section |] that in 
a QED plasma thermal corrections induce a non-Lorentz invariant 0(F 4 ) term which 
contributes to the photon splitting process. Compared to the 0(F 6 ) vacuum contribu- 
tion this thermal correction can actually be quite large already in an e + e~-plasma for 
T/m < 0.2, as indicated in Fig. [I], at least if eB/m 2 ~ 0.2. For an e~p + -plasma with 
chemical potential fi e g = m e there is a possibility of relatively large thermal correction at 
very low temperature (see Fig. 0) but the perturbative calculation in powers of eB breaks 
down in this limit. In Fig. |2| the Compton and photon splitting absorption coefficients 
become comparable for physical parameters corresponding to an electron density of the 
order of 10 28 cm~ 3 . In a neutron star such an electron density n e would occur only below 
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the surface of the star (the surface electron density can be the order of 10 cm ), where 



photon splitting processes are not likely to be of any importance anyhow [K|. In studies 
of one-temperature accretion disks around black holes temperatures like T ~ m can be 
achieved ||18|| . In such models the magnetic field is, however, very small (B ~ 10 5 G) and 
photon splitting is therefore most likely an unimportant physical process. We have not 
taken dispersive effects into account. For soft gamma-ray repeaters, dispersive effects can 
be of importance for large electron (and/or positron) densities in that e.g. the rate for 
Compton scattering is modified JTT 
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